CO 



O 



Van Hove Exciton-Cageons and High-T c Superconductivity: VIIID 
Solitons and Nonlinear Dynamics 

R.S. Markiewicz 

Physics Department and Barnett Institute, Northeastern University, Boston, MA 02115 
Running Title: Solitons and Nonlinear Dynamics 

Keywords: anharmonicity, charge- density waves, electron-phonon coupling, structural phase transition, soli- 
tons 



5_i 

& ' The low-temperature orthorhombic (LTO) phase transition in La2- £C Sr :E Cu04 can be interpreted as a 

dynamic Jahn- Teller effect, in which the degenerate electronic states are associated with the large densities 
of states at the two van Hove singularities. The equations describing this phase are strongly nonlinear. This 
paper illustrates some consequences of the nonlinearity, by presenting a rich variety of exact nonlinear wave 
solutions for the model. 

Of particular interest are soliton lattice solutions: arrays of domain walls separating regions of local 
low-temperature tetragonal (LTT) symmetry. These arrays have a macroscopic average symmetry higher 
CNI \ than LTT. These lattices can display either orthorhombic ('orthons') or tetragonal ('tetrons') symmetry, and 

can serve as models for a microscopic description of the dynamic JT LTO and high-temperature tetragonal 
£^ | phases, respectively. 

m , , „ 

i 1. Introduction 

d ' The two-dimensional van Hove singularity (vHs) was initially introduced as a model[l] to study competi- 

tion between superconductivity and structural instability. Its possible relevance to high-T c superconductivity 
was pointed out by Hirsch and Scalapino[2], who showed that, because of the logarithmically diverging den- 
sity of states (dos), the form of the BCS equation is modified in a vHs superconductor, leading to much 
higher values of T c . When the high-T c superconductors La2- a; Ba 2 ;Cu04 (LBCO) and La2- :E Sr :E Cu04 (LSCO) 
were first discovered, it was noted that their Fermi surface falls close to a vHs, and it was proposed that 
the vHs might be responsible for both the high T c values, as well as the structural transition from the 
high-temperature tetragonal (HTT) phase to the low-temperature orthorhombic (LTO) phase[3]. 
' However, the LTO phase was found not to split the degeneracy of the two vHs[4], and moreover did not 

have a large effect on the electronic properties of LSCO, so it was concluded that the LTO transition had 
nothing to do with the vHs. Now the absence of a structural transition does not rule out a role for the vHs 
in promoting superconductivity. Indeed, the T c enhancement depends mainly on the dos peak - i.e., on the 
q = charge susceptibility, and not the inter- vHs q = 2fcp susceptibility [5]. However, Phillips has shown[6] 
that a simple vHs model cannot explain many details of, e.g., the doping dependence of T c . Most of these 
problems could be overcome if the vHs also produced short-range structural disorder [7], but, according to 
Phillips, "there is no experimental evidence which would indicate any relation between the major lattice 
instabilities" and the vHs[6]. 

In fact, such evidence is starting to accumulate. First, there is a second structural phase transition 
in LBCO, this time to a low-temperature tetragonal (LTT) phase[8], which does have a large impact on 
electronic properties and essentially destroys superconductivity. Moreover, the transition is optimized at 
a fixed concentration of holes[9], and has the correct symmetry to split the vHs degeneracy [10]. Secondly, 
a careful reexamination of the LTO-HTT phase boundary has showed that superconductivity is clearly 
associated with the LTO phasefll]. Finally, it was recently proposed[12-14] that the LTO phase transition 
is indeed driven by the vHs, but the LTO phase is dynamic, and the average macroscopic symmetry is not 
representative of the local microscopic state. 

This new model describes the structural phase transitions in terms of the interplay between the two 
electronic saddle points (vHs's) and the optical phonons associated with tilting modes of the oxygen oc- 
tahedra, with subsidiary roles played by shear strains and acoustic phonons. The LTT phase involves an 
essentially static band Jahn- Teller (JT) effect[13], in which the two vHs's play the role of the degenerate 
electronic states. As the temperature is raised, there are large fluctuation effects, and the LTO and HTT 
phases can be described as dynamic and disordered JT phases, respectively [14]. 
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The model Hamiltonian is an extension of the Labbe-Friedel model for the A15 compounds[15,16]. It 
contains significant nonlincarities; the purpose of the present manuscript is to explicitly display some of 
these nonlinear features. Analysis of these nonlinear waves should prove useful in understanding the short- 
range order present in LBCO above the LTT transition, and as fluctuations in other high-T c cuprates[17], 
and in understanding the role of anharmonicity in enhancing T c . The present paper introduces a concrete 
microscopic model which displays the properties postulated for the dynamic JT phase. 

Many of the anomalous features are shared by a larger class of materials, the ferroelastics[18]. In addi- 
tion to the high-T c cuprates, this class includes a number of related perovskite compounds and ferroelectrics, 
including SrTiC>3, and BaTiC>3, and the martensitic phases of the A15's. Ferroelasticity is in general sug- 
gestive of collective dynamic JT effects, with strong electron-phonon coupling. Hence, the present model 
can be applied to transitions within this larger class of materials. Indeed, the model allows significant 
electron-phonon coupling even for a filled band, if the Fermi level falls in a gap between two hybridized 
bands [13]. 

The paper is organized as follows. Section 2 shows that corner-sharing of octahedra greatly restricts 
the possible low-lying states of the lattice, and requires the introduction of defects (solitons) to expand 
the variety of configurations. A possible model for domain-wall solitons is introduced, and it is shown 
how this state can provide a microscopic picture for the dynamic JT state of LSCO. In Section 3, the 
combined electron-phonon equations of motion are introduced, and reduced to a form which is convenient 
for analyzing the nonlinear wave solutions, both for waves propagating along 100 (x or y axes) or 110 
(at 45°) directions. Section 4 illustrates the variety of possible wave solutions, both dynamic and static. 
The low energy waves include both phason-like modes (weak, periodic modulation of the axis of tilt) and 
soliton-likc solutions. These latter can be either isolated domain walls between LTT domains, or arrays 
of solitons, which display an average macroscopic symmetry which is greater than the true microscopic 
(LTT) symmetry. The macroscopic symmetry can be either orthorhombic (the corresponding waves are 
called 'orthons') or tetragonal ('tetrons'). In Section 5, these results are applied to LSCO. The material 
parameters are estimated, and the energies of the resulting solitons are calculated. Soliton array solutions 
are predicted to arise for 100-type waves, including static arrays. The energy of individual domain walls 
is calculated, and compared to microscopic calculations (Appendix II). The relation of the various types 
of soliton arrays (tetron, orthon, isolated domain walls) and the observed phase transitions (HTT — > LTO 
— > LTT) in LSCO/LBCO is discussed. In Section 6, it is suggested that the solitons may be charged, and 
hence related to domain wall phases postulated earlier in doped LSCO. Finally, Section 7 summarizes the 
conclusions of this work. 

2. Model of Solitons 

The essential feature of the dynamic JT effect is that the global symmetry is not indicative of the 
local symmetry. Such a situation is found in other materials, in particular in the several phase transitions of 
BaTi03[19]. In the high temperature phases of BaTiOs, the sample is composed of microscopic domains of the 
low temperature (distorted) phase, but with different polarization axes in each domain. As the temperature 
is raised, the axes point successively over two, four, or all eight possible orientations, so that the sample as a 
whole appears to have higher and higher symmetry. These domains have been directly observed via streaking 
in electron diffraction patterns[20], and similar streaking has been observed in high-T c compounds[21]. 

The tilting in these materials has often been described at the mean-field level in terms of a pseudospin 
or Potts model, replacing the allowed tilts with equivalent spin states. A closely related mean-field treatment 
was applied to LSCO in Ref. [14] . A shortcoming of this procedure is that it overlooks the problems associated 
with corner or edge sharing of octahedra[22]. This corner sharing leads to a long-range interaction between 
tilts, and requires the introduction of macroscopic defects (solitons as domain walls) to break up the resulting 
LTT-type order. 

Let us illustrate this effect in LSCO. Consider a single octahedron, tilted about the x-axis (as in the 
LTT ground state). Since neighboring octahedra are corner-sharing, its nearest neighbors along the y-axis 
must be counter-tilted about the x-axis (again, as in the LTT phase). By induction, the whole row must 
be similarly tilted, in the absence of some form of excitation - either an anti-tilt along x or a crossover to 
a tilt about y. Similarly, if the entire row is tilted about x, the next row along the x-axis cannot be tilted 
about the y-axis, without tipping an oxygen which is shared with the original row. Hence, the second row 
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would also be tilted about the x-axis. By induction, every row would be tilted about the x-axis. Hence, in 
the absence of defects, the ground state would have a quasi-LTT pattern. To lowest order, there is only a 
weak interaction between the tilt direction on adjacent rows, so the low energy physics would (in the absence 
of defects) be that of rows of x-tilted octahedra, with two possible states for each row (e.g., any particular 
octahedron could have either a + or - tilt about x, which would fix the tilts of all other octahedra on that 
row). Hence, this physics would be in the universality class of the one-dimensional Ising model. 

Such one-dimensional chains have very low entropy: for an TV x N array of atoms there are N 2 atoms, 
but only N chains, each of which has only two possible orientations, for a total of 2 N configurations. To 
increase the entropy significantly, it is necessary to introduce macroscopic defects into the system. For 
instance, it would be preferable to have domains in half of which the tilts lie along the y-axis, in the other 
half along the x-axis. These domains of LTT-like phase must be separated by domain walls (solitons). The 
purpose of the present manuscript is to use the nonlinear equations of the system to develop a model for 
these domain walls. If the wells are deep, the domain walls are likely to be narrow - perhaps a single cell 
in width. If the wells are shallow (or at high temperatures) the walls can be broad as the octahedral tilt 
changes gradually. The average cell size will be given by the balance between surface tension and entropy. 

Figure 1 suggests a possible (narrow wall) soliton model of the LTO phase. The circles in Fig. 1 
represent the apical oxygens above the CuC>2 plane, which would be centered immediately above the Cu in 
the undistorted HTT phase. The arrows illustrate the instantaneous local distortion of the apical oxygen's 
position in the tilted state. The figure illustrates two LTT-like domains, separated by a dynamic domain 
wall, with rotating tilt direction, for motion either along a (110) direction (Fig. la) or along a (100) direction 
(Fig. lb). In the (110) wave, the domains are composed of diagonal stripes, alternating in tilt along the x- 
and y-axes. The line of atoms along the domain boundary cannot, because of corner sharing, be in a pure 
LTT tilt, but will assume a (dynamic) average tilt, as illustrated by the arcs of circles. If the octahedra 
within each cell have pure LTT tilts, the structure in the figure would not directly produce any LTO strain, 
e xy . However, there is a pronounced asymmetry to the tilts along directions parallel and perpendicular to the 
domain boundaries. Within each row parallel to the boundary, the tilt is uniquely either up or down, so that 
the entire domain can readily assume the corrugated pattern characteristic of the LTO phase (i.e., strong 
coupling to a e xy strain. On the other hand, in the perpenducular direction, the up and down tilted rows 
are interchanged each time a soliton wall is crossed, making it impossible to generate a global corrigation in 
this direction. Moreover, the octahedra comprising the soliton already lead to a net shear strain in the easy 
direction. 

Thus, the proposed configuration (1) has domains of local LTT symmetry, (2) arranged in such a 
pattern as to readily couple to a corrugating shear along a preferred easy axis, to give a macroscopic LTO 
symmetry. Finally, (3) the domain walls already provide the appropriate shear strain, so that e xy should be 
approximately proportional to the density of solitons. 

It is very often convenient, in dealing with structural phase transitions, to think in terms of pseudospins, 
since the magnetic analogs are generally much better developed. Within such a pseudospin model, the solitons 
have a very natural interpretation: the two kinds of domain correspond to S z > and S z < 0, while the 
soliton is simply a Neel wall, with S x ^ inside the wall. 

3. Real-Space Hamiltonian: Intracell and Intercell Coupling 

3a. Tilt Equation of Motion 

In the present paper, the Hamiltonian will be a simplified form of that introduced in Ref. [14] (Eq. 
25 of that paper). The coupling to strains and umklapp electron sacttering were found to be relatively 
unimportant, and hence will be neglected. Furthermore, it will be assumed that the octahedra are all tilted, 
but in random directions, so the pseudorotation operators have the form R\ — Rcos<j>, R2 = Rsinifi, with R 
constant. The Hamiltonian then becomes 

I 

+ Fo E[( ^(0-^(r+^) )4 + { R 2 (l)-R,(Uvax) )4] + T_l [m + ^ + T, Rim ^ (lfl) 
v 



3 



^ -Ri^-Ei^ + ^ay) 2 i?2(0 - J? 2 (f+ 7?ax) 2 
W ep ; = -a_| 2^K 2 ) ~( 2 ' ^ ' 

n 

where n is summed over ±1. It should be noted that R4 does not have the dimensions of a lcngth[13]; rather 

where m is the mass of an oxygen atom and Z is the displacement of the planar oxygen above/below the 
Cu0 2 plane caused by the pseudorotation R t . This means that the parameter M is dimensionless, 

M=±(l + rj 2 ap )~3M, 

where r\ ap ~ 1.27 is the ratio of the Cu-apical O distance to the Cu-planar O distance. 

In Eq. lb, H ep i is the electron-phonon coupling term in the Zth cell, and the absolute value arises from 
the JT effect [14]. This absolute value sign leads to unnecessary complications in the analysis, since the model 
is only approximate to begin with, so a further approximation is introduced. This approximation essentially 
amounts to Fourier analyzing the electron-phonon term in cos{Aruf)), and retaining only the lowest order 
terms. This procedure is carried out in Appendix I. 

The formulation of Eq. 1 is similar to calculations introduced to describe the anomalous polarizability 
of 2 ~ in ferroelectrics[23,24], and the resulting nonlinear waves are also similar to those found in fcrro- 
electrics[25]. The equation of motion corresponding to Eq. 1 then becomes: 



"52^(0 + "7T \^[cos<j>(l)(siruf>(l) — sin4>(l + rjax)) 3 — siruf>(l)(cos(j>(l) — cos<j>{l + r/ay)) 3 ] 
R 8 

v 

+ (2T a Q - T 2 )sin<}>{l)cosci){l)[sin<j> 2 {l) - cos<f> 2 (T)} + w — % = 0, (2) 

-K d(j)(l) 

where the last term is given by Eq. A4 in Appendix I. 
3b. (110) Waves 

In studying the nonlinear wave solutions of Eq. 2, it is convenient to look for one-dimensional solutions. 
The simplest such solutions propagate along the x or y axes or at 45°. It is assumed that the background 
state is of HTT symmetry (i.e., sin<j) — cos(f> — except for the modulation produced by the wave), and that 
the tilt has the same magnitude for all atoms in any plane perpendicular to the propagation direction, but 
with neighboring cages tilting in the opposite sense. As suggested by Fig. 1, the most interesting solutions 
are those propagating at 45°, (110) waves, which include the LTO-like domain walls, and this case will be 
analyzed in the present section. 

Figure 1 can be used to understand the reduction of the problem to one dimension. For convenience, 
lines of octahedra parallel to the propagation direction (e.g., A and B in Fig. 1) will be referred to as 
columns, and lines perpendicular to this direction, forming wave fronts, as rows. Thus, it can be seen in Fig. 
1 that all atoms fall into one of two classes of columns, A or B, and likewise all rows, and that all of the 
nearest neighbors of an A octahedron are B octahedra. In a one-dimensional, 45° wave, it will be assumed 
that all the octahedra within a given row have exactly the same value of (j>, so that the <j>(l) may be rediced 
to the one dimensional <f> n , where n denotes the successive row in the propagation direction. Because of 
the antiferrodistortive nature of the transition, the tilts of adjacent cells tend to be opposite in phase. It is 
convenient to explicitly account for this factor, by defining 

cos(j)(l) = {-l) n cos(j) n . (3) 

The problem is thus effectively reduced to one-dimensional. For the 45° waves, Eq. 2a takes the form 



h ^ (2r a r 2 ) 2 oh 

-jg&n + y^}cos<j) n {sin(j) n - sin(j} n+ri ) 3 - sin<j> n {cos(]) n - cos<j) n+11 ) 3 ] ^ — —sinA<f> n + g * pn = 0. 

(4) 
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Whereas anharmonic phonons have often been studied as a paradigm for soliton physics, most models assume 
a diagonal anharmonicity. In contrast, the present model has significant nonlinearity also in the intercell 
coupling. 

Equation 4 has travelling wave solutions of the form: 

4> n+1 {t + a'/v) = <j> n {t), (5) 

where a' = a/y/2 is the effective repeat distance, a is the Cu-Cu distance and v is the phase velocity. The 
continuum limit amounts to expanding 

sincj) n±1 (t) = sin<t>„(t) ± TCOs<j> n {t)j> n (t) + ^T 2 (coscf) n {t)4> n {t) - sin<j> n (t)<f>l(i)), (6) 

where r = a'/v, with a similar expansion for cos(f>. Substituting Eq. 6 into Eq. 5 reduces that equation to 
a single-site form: 

m (2r + 2r a r i 3T 2 r 

2_ V — 2__2 ^ S m40+^^[0 2 sin4(/.+2^sm 2 20]+r ep [4sm40+T 2 (^cos40-0 2 sm40)] = 0, (7a) 

R 4 8 

where terms of higher order than r 2 have been neglected, and the index n has been omitted. The final term 
of Eq. 7 follows from Eq. A6 of Appendix I, with 



and 

in which case Eq. 7a becomes 



Equation 7 may be rewritten in dimensionless form by introducing 

p , = (i6r eP + r 2 -jr - 2 rg)^ ^ 

3rW = ± 

7 32M 4w 2 ' 1 ' 

a = 40, (8d) 
+ 7« 2 )sma 

a = — i : 1 ( 9 ) 

1 + A'ycosa 

with /3 = /3'/(l — 47'), and 7 = (S 1 — — 47'). The quantities 7' and <5' are intrinsically positive. In 

principle, (3' (Eq. 8a) could be either positive or negative, but will here be taken positive, to stabilize an 
LTT-like phase (potential minima at (j> = 0, 7r/2, etc.). However, (3 in Eq. 9 will be negative when 47' > 1 
(v < vi), while 7 may have either sign. When the JT coupling effects dominate over the anharmonicity, then 
5' > 7', so 7 will be positive for fast waves. The velocity v\ is a critical velocity, which separates fast waves 
(v > vi) from slow waves (v < vi). 

3c. (100) Waves 

The analysis for (100) waves is very similar to the above analysis for (110) waves. In this Section, a 
subscript will be added to various symbols to distinguish 100 from 110 waves. Since both waves satisfy the 
same equation, the symbol will be omitted in subsequent sections, unless it is necessary to distinguish the 
two types of wave (e.g., in Section 5). Now the wave propagates along the x-direction, with rows along y. 
Hence, nearest neighbors along the rows are exactly out of phase, 

R 1 (T±ay) = -R 1 (T). 
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By a calculation similar to that given above, the (100) wave equation has the same form as Eq. 9, with 
the following changes. First, ti o = a/v, since a is the distance between rows. Next, since only half of the 
neighbors now contribute to the dispersion, 

Tioo _ 7no 



t 2 9t 2 ' 

'100 ZT 110 

with a similar relation for S' 100 . However, since a' 2 = a 2 /2, r 2 00 = 2t 2 10 , so 7^ 00 = 7 110 and S[ 00 = S' 11Q . 
Thus, these two corrections compensate. However, there is one last correction, an additional term in H epn , 
Eq. A9 in Appendix I. This changes the denominator in the expressions for both (3 and 7: 

(l-4 7 ')-(l-4(7 / + <5 / ))- (10) 

Hence, (100) waves have virtually the same dispersion as (110) waves. Both waves satisfy Eq. 9, with 
the only change being the substitution, Eq. 10, for the (100) waves. In Section 5, the significance of this 
difference will become apparent. Section 4 analyzes possible solutions of Eq. 9, and hence is applicable to 
both types of waves. 

4. Nonlinear Wave Solutions 

4a. Phase Portraits 

A first integral of Eq. 9 can readily be found by setting y = a, 

dy . , 

a = y = —a = y y, 
da 



so 



' yo 

or 



f y ydy f a sinada 

Jyo IV 2 ~ P Jao 1 + 4 7cosa ' 



V2 - P - l 1 + ^ C0Sa °^F(a), (11) 



y 2 - P V 1 + ijcosa 

where «o and yo arc initial values of a and y respectively, and $ — (3/-f. Equation 11 may be reduced to 
quadrature: 

t= / = . (12) 

J °o sj F{a)(yl - (3) + (3 

This equation must be solved numerically. Equation 11 depends on two parameters [3 and 7. In constructing 
phase portraits, only 7 > need be considered, since 7 < is equivalent to 7 > 0, but with a — > a + n. The 
nature of the phase portraits does change, however, when 7 = 7 C = 1/4. Hence, there are four classes of 
phase portrait, corresponding to (3 >, < 0, and 7 >, < "f c . For LSCO, the fast wave solutions are expected 
to correspond to (3, 7 > 0, and slow waves to /3, 7 < 0. The corresponding phase portraits (a vs a) are 
illustrated in Figs. 2 (/3,7 > 0), 3 {(3 < 0, 7 < 0), and 4 ((3 < 0, 7 > 0). The lines with filled circles in these 
figures are separatrices, dividing regions of qualitatively different kinds of wave behavior. The following 
subsections will describe the various kinds of waves which can be generated. 

The situation corresponding to Fig. 4 does not appear to arise in LSCO. However, it could occur 
under special circumstances, and since the resulting waves are strikingly different (particularly for 7 > 7 C ), 
solutions for this case will be discussed in Section 4d. 



4b. Fast Wave Solutions 

For (3 > and < 7 < j c , the typical phase space profile (a vs a) is shown in Fig. 2a. This diagram 
periodically repeats outside the range illustrated, a — > a + 2nir, n integer. The overall nature of the phase 
portraits can be understood by noting that, for 7 — > 0, Eq. 9 reduces to the equation of a pendulum. Hence, 
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there will be low energy solutions corresponding to oscillations about a single potential minimum, and higher 
energy solutions in which the octahedron hops between adjacent minima (the 'pendulum' undergoes complete 
2ir rotations). The finite values for 7, due to intercell coupling, lead to all the complications discussed below. 

In Fig. 2a, there are two separatrices, the nature of which can best be understood by looking at a series 
of typical orbits, Fig. 5. These orbits all correspond to the parameters /? = 7 = 0.2, a = 0, with different 
choices of yo = a(t = 0). This corresponds to moving up the left hand axis of Fig. 2a. Here and below, the 
waveforms are found by numerical integration of Eq. 12. For the parameters of Fig. 5, the first separatrix 
occurs at y — yi, with 0.816 < y\ < 0.817. For y < yi, Fig. 2a shows that there are a values which are not 
sampled (the phase space plot passes through a = 0). The resulting orbit is consequently a bounded periodic 
orbit. For yo > y c , all values of a are sampled, and, since a never vanishes, a increases monotonically. 

An analytic formula for y\ can readily be found from Eq. 11. From Fig. 2a, the separatrix has initial 
conditions a = 0, y = j/i, and passes through the point a = n, y = 0. Then Eq. 11 may be rearranged to 



When a = n — a' with a' « 1, the integrand is <~ 1/a'. Thus, the integral diverges logarithmically: the 
separatrix is a 'single kink' soliton, varying once from — it to +ir as t varies from —00 to +00. 



The second separatrix corresponds to yo = 2/2 = y $■ From Eq. 9, it can be seen that a — when 
yo = ?/2- Hence, a = y 2 remains constant, and a increases linearly with time. For y > y 2 , a again increases 
monotonically with time. However, now a is larger near a = (2n + l)ir, than near 0, 2mr, so a tends to 
linger near the potential minima (lower slope near 2nn). In contrast, for y < 1, the orbit tends to avoid the 
potential minima. The extreme case is for the separatrix y = yi, where the orbit is a single kink, Eq. 13. 
These high-d solutions are not of much interest, since they also correspond to high energies; the low energy 
solutions are those which spend most of their time near the potential minima. 

Figure 6 shows a number of different periodic orbits, for y < y\. The low energy waves are the small 
amplitude, nearly sinusoidal waves which remain close to a = 0. These correspond to waves in an LTT-type 
phase: the octahedra are all tilted in nearly the same x-direction ((f) = a/ A ~ 0, but the tilt axis wanders back 
and forth as the wave progresses. As the amplitude increases, the waves become increasingly nonsinusoidal, 
tending to spend more time near the potential maxima a = ±7r. These are not of importance, being high 
energy waves, but qualitatively similar slow waves are of great interest, and are discussed further in Section 
4c. 

4c. Slow Waves: Soliton Lattices in the 'LTO Phase' 

When v becomes less than v\, the signs of both [3 and 7 change. This sign change is completely 
compensated for by shifting a — > a + 7r, as can be seen by comparing Figs. 2 and 4. Thus, the periodic orbits 
of Fig. 6 can be interpreted as fast wave solutions, using the left hand axis, or cquivalently as slow wave 
solutions (with (3, 7 of same magnitude but opposite sign), using the right hand axis. This simple switch 
has a profound effect on the nature of the low energy solutions. Now the small amplitude waves are high 
energy waves, corresponding to small oscillations about the unstable energy maximum <f) = tt/4. 

On the other hand, the large amplitude waves arc now low energy solutions, because they spend most 
of their time near potential minima at <j> = or ir/2. These are soliton lattice solutions, and a particular 
solution is displayed in Fig. 7. Corresponding to the periodic solutions for y < y\, there are unbounded 
solutions for y > yi- The most interesting of these are ones, such as illustrated in Fig. 8, which again form 
soliton lattices. 

Since these are travelling wave solutions, Eq. 5, the spatial dependence can be restored by substituting 
t —* t — x/v, where in the lattice case, x — na', with n an integer. Hence, Fig. 7 can be thought of as an 
illustration of the spatial variation of at a fixed instant of time. Finally, in the limit of a static solution 




(13a) 



which yields y\ — 0.8165 for = 7 = 0.2. The time to travel from a = to a = tt is 




(136) 
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v — ► 0, Fig. 7 would be the spatial profile of a static domain array, transforming from domains of x-tiltcd 
LTT phase ((f) = 0) to domains of y-tilted LTT phase (<f> — n/2), precisely as in Fig. 1. 

The bounded orbit solutions (Figure 7) may be called 'orthons', in that a static array of orthons would 
produce a phase of average macroscopic orthorhombic symmetry - the material keeps switching back and 
forth between an x-tilted domain and a y-tilted domain. The unbounded solutions (Figure 8) are 'tetrons' 
in that they switch periodically between domains of tilt +x, +y, —x, —y, +x, to produce an average 
tetragonal symmetry 

The distinction between orthons and tetrons can be more clearly understood by analyzing the associated 
strains. From Ref. [13,14], it is known that the orthorhombic shear strain e xy is proportional to the product 
R1R2, whereas the LTT phase (in a single layer) is associated with e_ = en — 622 oc R\ — R 2 . Figure 7b, c 
(8b, c) plot the values of cos2<j) and sin2<p, which should be proportional to e_ and ei2, respectively, in the 
orthon (tetron) wave. Within the individual domains, there is a large value of e_. However, the sign of 
e_ alternates in successive domains, so that the macroscopic average value of e_ vanishes. In contrast, the 
octahedral shear strain en is non-zero only within the domain walls. However, for an orthon, this strain has 
the same sign in every domain wall, so that there is a net macroscopic shear < ei2 >, proportional to the 
domain wall density. For a tetron, the shear strain is equally likely to be positive or negative. Thus, despite 
the presence of many LTT- type domains, the tetron has on average no macroscopic shear strains, leading to 
a macroscopic HTT symmetry. 

4d. 7 > 7 C 

When 7 > 7c, Figs. 2b, 3b, 4b, the slope a passes through a critical value at special angles a c such 
that cosa c = I/47 (Eq. 11). For (3 > 0, a diverges when yo > (3 (case A) and a —* when yo < f3 (case B), 
leading to strikingly different behavior. In case A, the resulting soliton lattices are qualitatively similar to 
those found for 7 < j c , except that a jumps discontinuously from a c to 27r — a c . Similar behavior arises when 
(3 < (case C), as illustrated in Fig. 9, for a tetron lattice. (For case C, orthon arrays with a discontinuities 
are also possible.) This discontinuous jump means that the interface is abrupt, and its properties must be 
calculated on a microscopic model, and not in the continuum limit. 

For case B, the solutions are again periodic oscillations, similar to those in Fig. 6. However, the lines 
a = a c , 2n — a c divide these oscillations into those centered on 2nn and those centered on (2n + l)w. 

Soliton arrays with jump discontinuities are an intriguing possibility, but are found only in case C. 
While such a situation does not appear to arise in LSCO for 100 or 110 waves, it might still occur in other 
directions or in materials with different parameters. Thus, if anharmonic effects were stronger, it would be 
possible to have [3' > 0, but 5' < 7', so j3 < 0. 

4e. Static Solutions 

Some care must be exercised in examining the static limit of Eq. 9. The static solutions of Eq. 9 
correspond to v — > 0, or equivalently t — > 00; however, the time derivatives are really derivatives with 
respect to t — x/v. Hence, changing to spatial derivatives requires multiplying a or a 2 by v 2 before taking 
the limit. Hence, the static limit of Eq. 9 still has the form of Eq. 9, but with the superscript dots 
now signifying spatial derivatives, with [3 —f3'/A"f'v 2 and 7 — > —(5' — j')/4r/'. That is, both (3 and 7 
are expected to be negative in this limit. In this case, the phase profiles are identical to the slow waves, 
discussed in Section 4c. 

4f. Energy of a Soliton 

While the above analysis shows what kinds of nonlinear waves are possible, no account was taken as 
to the energy of the wave, to see which waves are probable. The low energy waves will be those which are 
close in energy to the ground state of the system - i.e., to the LTT phase with <f> a multiple of n/2. These 
low-energy waves can be of two forms. First, there are ordinary phasons - nearly sinusoidal excitations of 
small amplitude away from the LTT phase, as in Fig. 6. 

More interesting, however, are the soliton lattice solutions, Figs. 7-9. These will be low energy states 
when the steps are flat topped and lie close to an LTT value, <fi — nw/2, for integral n. Orthon solutions of 
this form can be generated from Eq. 11 by chosing y — and a small. As a gets progressively smaller, 
the steps become flatter and the periodicity increases. In the limit a — > 0, the solution approaches the 



8 



separatrix, which is an isolated soliton. In the opposite limit, choosing a = and yo small generates low 
energy tetron arrays. 

The energy per unit cell associated with the wave can be found from Eq. 1. Applying the same 
transformations as to the equation of motion, the energy is found to be 

MR? 

E = -^—4> 2 - {Pe - 7e^ 2 )cos4(/) + 5 e '(j)sinA<j), (14) 



with 



D, = ( ^ + r i -2rg-2r„ )fl , | (15o) 

S. = F'-KW . (15c, 
M = M — 3r »f T \ (15d) 



The static limit of Eq. 14 follows from the substitution 

d d ttfi\ 

T dt^~Te (16) 

with £ = x/a! (for (110) waves). Denoting the spatial derivative on the right side of Eq. 16 by a prime, Eq. 
14 becomes 

E = E (-l<j)' 2 - e - %<j) ,2 )cos^ + 6 e <j)"sin4<j)), (17) 

with 

E = T R\ (18a) 
K = S ^ + ^~ 2 , (186) 

7e = (18c) 

Se = (18d) 

and ?7 e p = r ep /r , r\i = (r 2 — 2rg)/r . The static limit of the parameters of Eq. 9 can also be expressed in 
terms of the 77 parameters. However, these differ for 110 and 100 waves: 

g 

/?no = -g(16»7e P + »j2-2), (19a) 
7iio = -^(8% P -3). (196) 

*»~i(£^§)- < 1M > 

Note that /3 remains the same for both kinds of waves. 

Figure 10 illustrates the energy per unit cell, Eq. 17, for the oscillatory sloutions of Fig. 6, assuming 
i] ep = 2, r/2 = 0. [N.B.: This is only done for illustrative purposes. Strictly speaking, the assumed r\ values 
are not consistent with the parameters of Fig. 6; this will be corrected in Section 5, after the parameters 
relevant to LSCO have been estimated.] For simplicity, the energy associated with only one half-cycle of 
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oscillation (i.e., a single soliton) is explicitly displayed. In this case, the soliton has a simple lineshape, the 



energy being dominated by the term (3 e (Eq. 18b). The linewidth is ~ 5a/y /?, and the profile quickly 
becomes independent of soliton spacing - i.e., there is very little soliton-soliton interaction in the present 
model. 

5. Applications to LSCO 

Numerical estimates for the above parameters follow from the calculations of Ref. [13]. These calcula- 
tions may underestimate the parameters, since it was assumed that the HTT^LTO transition represents the 
onset of tilting, R^ 0, whereas it is more probable that R is finite in the HTT phase, but with < coscj) >=0. 
Furthermore, the value of R is taken from the measured tilt, whereas this macroscopic average value may 
only be about half as large as the instantaneous microscopic value [14]. With these caveats, the best estimates 
are from parameter sets 10 and 11 of Table II of Ref. [13]. The resulting parameter values are summarized 
in Table I. The Table also includes values for the crossover velocity, v\, where 



Note the striking difference of 7 values for the two types of waves. For 110 waves, 7 is large and negative, 
dominated by the electron-phonon coupling. However, for 100 waves, \^\ is always restricted to be less than 
7 C = 1/4. Hence, the slow wave and static solutions for 100 waves are described by the phase profiles of Fig. 
3a, and include orthons, tetrons, and isolated kinks. By contrast, the 110 waves correspond to Fig. 3b, so 
the low energy solutions are phasons: weak periodic modulations of the tilt axis about some LTT state. 

Figure 11 illustrates the energy of a soliton wall, using Eqs. 18 and 19 to estimate all parameters, and 
using typical values rj ep = 8, r\i = 0, and Eq = h.bmeV. Shown are both the total energy (solid line) and 
the various components (corresponding to Eqs. 18a-d). The shape is more complicated than found in Fig. 
10, because the term in % (Eq. 18c) becomes comparable to that in $ e when \$\ is large. 

By integrating the energy profile, the total interface energy E w (plotted as energy per column width) 
can be found, and is plotted in Fig. 12 as a function of the spacing between successive domain walls, d w . In 
the present theory, the tilt-tilt interaction is short ranged, and this is reflected in Fig. 12 by the fact that 
there is virtually no soliton-soliton interaction until the solitons are less than two cells apart. Below this 
separation, the tetrons have lower energy than the orthons. However, it is unlikely that a continuum theory 
should be trusted at such small length scales - note that the width of the domain wall is less than a single 
cell across. The limiting energy of a single soliton at large separation is found to be nearly linear in i] ep : 



There are a number of additional factors which will give rise to longer range soliton interaction. For 
instance, the strains associated with the tilting give rise to direct coupling of more distant cclls[14]. Moreover, 
the electronic energy lowering is produced by the splitting of the vHs degeneracy within the LTT-like domains. 
If the domain walls approach one another too closely, the domains will cease to have a well-defined LTT-like 
character, and the vHs splitting will decrease. Thie is equivalent to a repulsive interaction between domain 
walls. 

Using the present equations, it is possible to calculate the energy associated with a static domain wall 
exactly, for a given microscopic configuration, Fig. 13. Such calculations are presented in Appendix II. 
These calculations are still approximate, in that the simple distortions assumed do not display the correct 
corner sharing behavior. The resulting energy is very similar to the calculations of this section. The wall 
energy is found to satisfy an equation of the form of Eq. 20, but with different numerical coefficients, which 
depend explicitly on the assumed wall thickness. The minimum wall energy corresponds to an abrupt wall, 
as in Fig. 13a. 

It is tempting to relate the various nonlinear waves to the various phases found in LSCO and related 
compounds - the isolated solitons separating large grains of LTT phase, the orthons associated with the 
LTO phase, and the tetrons with the HTT phase. If this correspondence is valid, then what distinguishes 
the various phases, and in particular, what constitutes a phase transition? Presumably, this is accomplished 





3E a 2 
16MR 2 ' 



E w ~^-( Vep -1.3). 



(20) 
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by coupling to macroscopic strains: e.g., in the LTO phase, the macroscopic average shear strain, e xy . As 
shown in Ref. [14], inclusion of strain coupling makes the effective tilt-tilt interaction longer ranged, and 
hence could stabilize an orthon or tetron phase. [This same model could be applied to BaTiOs, in which 
case, the long-range coupling would be via the ferroelectric polarization.] 

6. Charged Solitons? 

In the above sections, it has been shown that the LTO phase of LSCO can be interpreted as a dynamic 
JT phase, or equivalently as a soliton lattice with domains of LTT phase separated by domain walls of a 
different phase, closer to the usual LTO phase. This nanoscale phase heterogeneity is very reminiscent of 
another grain boundary phase postulated[7,26] to exist in this system - a phase separation of holes associated 
with doping LSCO away from the optimal hole concentration at which the vHs coincides with the Fermi 
level. 

In this earlier model, the phase separation arises because the electron-phonon coupling is strongest 
exactly at the vHs, leading to a minimum in the free energy at precisely that concentration. (There is a 
close relation to Hume-Rothery phases in alloys.) The size of the domains is restricted, because the hole 
concentration is different in the two phases - i.e., the domains are charged, and the domain size is limited by 
Coulomb effects. The lowest energy phase turns out to be a grain boundary phase, with the LTT domains 
forming one phase, and the domain walls constituting a second, microscopic LTO-like phase. Near the vHs, 
it is possible to have a pure LTT phase. As the material is doped away from the vHs toward half-filling, the 
LTO phase becomes more stable, and the size of the LTT domains decreases. If, as found above (Fig. 7), 
the orthorhombic shear strain is taken as proportional to the domain wall density, this is exactly what is 
found experimentally in LSCO - that ei2 increases as the material is doped away from the vHs (x decreased 
below - 0.15). 

There are a number of similarities between the doping-induced grain boundary phase and the soliton 
lattice. Indeed, it is possible that the solitons are charged. This could arise if the strains in the domain wall, 
associated with corner sharing, lead to a modification of the shape of the octahedron - i.e., a change of the 
Cu to apical O bond length. A reduction of this bond length would reduce the splitting between the Cu 
d 3z 2_ r 2 and d x 2_ y 2 levels, which can lead to a different degree of charge transfer into the Cu-0 2 planes, as 
well as to a direct change in the shape of the Fermi surface. These changes will have similar effects to the 
charge bunching postulated in the earlier work[26], in that the density of planar Cu holes will be different in 
the domains and in the walls. It is of interest to note that a number of experiments have detected structural 
anomalies in the cuprates, often at temperatures very close to T c [17], which are related to tilting of the 
octahedra, and generally include a change of the Cu to apical O bond length. 

7. Conclusions and Discussion 

The present paper is the fourth in a series of closely related works which have introduced the idea that 
the LTO phase of LSCO is a dynamic JT phase. It is perhaps appropriate to briefly summarize what has 
been accomplished to date, and what remains to be done. The formal model was introduced in B[13], and 
static JT solutions were displayed, emphasizing the first appearance of a macroscopic tilt, R ^ 0. The 
dynamics of the model, wherein the tilt axis could bob about, assuming one of four preferred orientations, is 
discussed in C[14] at the mean-field level, and in D, the present manuscript, in terms of a microscopic picture 
of nonlinear wave solutions. The first paper, A[12], introduced the concept of the dynamic vHs-JT phase, 
and provided a simplified model calculation which reproduced the general doping dependence of the phase 
diagram. In particular, it answered the vexing question of why, if superconductivity can only be found in 
the LTO phase[ll], the highest superconducting transition temperatures occur in the doping range in which 
the LTO phase is least stable. 

Prior to these papers, the vHs model seemed to provide an interesting framework for interpreting many 
of the phenomena associated with high-T c superconductivity [7, 27], but also predicted a strong competition 
with structural instabilities (charge density wave phenomena) of a form which could not be clearly identified 
with the major phase transitions of the high-T c cuprates. While the discovery of the LTT phase provided a 
first link between structural instability and the vHs, the present model clearly demonstrates for the first time 
a close relation between the vHs and the LTO phase, the dominant phase in which high-T c superconductivity 
arises. Moreover, the model paints this LTO phase as a highly unusual, dynamic and nonlinear phase, of just 
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the form wherein one might expect to observe the anomalous normal-state and superconducting properties 
actually found in the cuprates. The close similarity between the present model and the earlier picture of 
short-range structural disorder[7,26] should be noted. 

While the present equations display a rich variety of nonlinear wave solutions, it should be kept in mind 
that these equations are already approximations, and that the more exact equations[13] should have even 
stronger nonlinearities. In particular, it was assumed that the tilting of the individual octahedral cages sets 
in at considerably higher temperatures, so that the average tilt R can be considered constant. Moreover, 
while the term H ep in Eq. 1 accounts for the role of the electronic excitations in modifying the structural 
transition, the back-reaction of this structural change on the electronic properties has not been discussed. 
In particular, what is the effect of the domain walls on the splitting of the vHs? Indeed, in a dynamic JT 
phase, it is unclear what the 'band structure' means, when the Brillouin zone itself may be a local function 
of time and space (c.f. the spatial variation of the orthorhombic splitting in Figs. 7 and 8). 

It should be noted that these complications are 'interesting' in that they are likely to lead to a rich 
variety of anomalous behavior, both in the superconducting and normal states, comparable to that observed 
experimentally in the cuprates. 

Acknowledgments: I would like to thank J. Jose, A. Widom, and F.Y. Wu for stimulating conversa- 
tions. Publication 546 from the Barnctt Institute. 



Appendix I: Approximate Electron-Phonon Coupling Term 

In this Appendix, an approximate form for the electron-phonon coupling term, Eq. lb, will be derived. 
Following the discussion of Section 3, this term can be written 

ot e R? 

H epn = — \^[(cos<j) n - cos<j) n+v ) 2 - (sin(f> n - sin(t) n+n )) 2 ]\. (Al) 

v 

Expanding sincf) n+71 , cos(f> n+n as in Eq. 6, this becomes 

a e _R 2 

H epn = —\Acos2cj) - T 2 (2<j)sin2(f) + 3<i> 2 cos2<j))\, (A2a) 

H epn = ^—(4 - 3t 2 2 ) v /1 + M 2 |cos(2 ( / ) + u)\, (A2b) 



where the index n on <fi is omitted, /i = 2r 2 </>/(4 — 3r 2 ), and sinv = fi/\/l + [x 2 - Figure 14a illustrates 
H epn in the case of perfect intercell coupling, cos0„±i = —cos(j) n (0 = = 0). Deviations from this coupling 
reduce the amplitude of H epn and shift the phase of cos2cj). 

The absolute value signs in Eq. A2 lead to the cusps in Fig. 14a, and can complicate the analysis. Since 
H epn is an effective potential for the tilt phonons, rounding the cusps at the potential maxima should have 
negligible effect on the problem. Hence, the Fourier expansion may be utilized: 

2 00 2(— F)™^ 1 
|cos20| = -(1 + V / cosAn4>). (A3) 

n—l 

Keeping only the constant and n = 1 terms provides a reasonable first approximation (dashed line in Fig. 
14a). 

What enters into the equations of motion is not H epn but its derivative with respect to (fi n , which may 
be written 

dH a. e R 2 

— = (±)— — [4sm20 + T 2 ((j)cos24> - 2 sin2<f>)], (AAa) 

0(p n z 

= (±)?^(4-T 2 l j ) 2 )^TT^sin(20 + v'), (A4b) 

where, to order r 2 , fi' = [i and v' = v, and the term (±) comes from the absolute value signs in Eq. Al: the 
minus sign applies whenever the quantity within ||'s in Eq. Al is negative. The resulting quantity, plotted 
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in Fig. 14b, has discontinuous jumps whenever the cos 2 (/>-term changes sign, Fig. 14a. Because of these 
discontinuities, the Fourier series requires more terms for an adequate representation: 

8 °° n(-l)" +1 

(±)sin2(f) = - ^ 2 ! sinAmp. (A5) 

n=l 

Figure 14b shows the Fourier series for one term (dashed line) and for two (dotted line). The one term series 
is nonetheless qualitatively correct, particularly in the region of the potential minima, so this approximation 
is adopted for the main text. In the present case, this amounts to replacing 

g 

(±)sin(2(f> + v) -> —sin(Ancj) + v), (A6) 

SO 

dH 4:6t e R 2 

-— ^ = — - — [AsinAcf) + t 2 {4>cosA<t> - <j) 2 sinA(j>)\ . (A7) 
o<Pn 3ir 

There is one last complication. The terms i? ep ,n±i also depend on <f> n , so their derivatives must be 
added to Eq. A7 to give the total potential contribution. By reasoning similar to the above, it can be shown 
that 



(A8) 



except possibly when cos2(f> n is changing sign. In this case, it is possible that, e.g., dH epn+ i/d<j> n will change 
sign before dH epn /d<f) n (Fig. 14b). In this case, its contribution would tend to cancel rather than add. This 
effect will tend to smooth out the discontinuity in Fig. 14b. Furthermore, the correction will be small in the 
smoothed form of the restoring force (dashed line in Fig. 14b). Hence, this complication is neglected, and 
Eq. A8 is assumed to hold throughout. This amounts to doubling the contribution of H epn , accounting for 
a factor of 2 in Eq. 2 and subsequent equations. 

Finally, the above analysis can be repeated for (100) waves, as in Section 3c. In this case, Eq. A7 
becomes 

= -^—[4sin4 ( p+—U(cos4<l>+l) - 2 szn40)l. (A9) 
oq> n 3tt 2 

Note that this differs from Eq. A7, in that the r 2 terms are half as large, and there is an extra term 
multiplying <fr. 

Appendix II: Microscopic Domain Wall Energy Calculations 

For narrow domain walls (only a few unit cells thick), it is possible to directly calculate the domain wall 
energy from Eq. 1, for any assumed pattern of tilt angles. In this Appendix, the calculation is presented for 
the three configurations illustrated in Fig. 13. In this figure, the dashed lines enclose the regions of canted 
cells - everything outside the dashed lines is assumed to be in perfect LTT order. (For the abrupt wall of 
Fig. 13a, there are no cells 'within' the wall.) The domain wall energy is defined as the excess energy over 
the ground state energy within a single domain. For an LTT- type domain, e.g., cos0 = everywhere, this 
minimum energy is, for a single cell, 

E mm = (^±^)R*~2^R 2 . (Bl) 

(Recall that the terms involving Tg and T2 are on-site terms, while the terms in Tq and a e _ involve sums 
over nearest neighbors.) 

For the abrupt wall of Fig. 13a, all atoms remain tilted about the x or y axes, so there is no wall energy 
contribution from the terms involving T2 or Tg. Here, y is defined as parallel to the domain wall and x 
perpendicular to it, with the positive x direction pointing to the right of the figure. The terms in To and a e _ 
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are affected by the wall, since the cells nearest the walls have neighbors in a less than optimal configuration. 
For instance, for the top left atom in the figure, the excess energy due to r is 

A J E r = E^(2xl + 0+(i) 4 -2). 

Here the first term in parentheses is due to the term involving nearest neighbors along y, the next to the 
neighbor along -x, the third along +x (i.e., across the domain wall); the final term is the subtraction of the 
energy for the undistorted state (absence of the domain wall) . For the top right atom, the analogous term is 

A£ r = M!(2xO + l + (i) 4 -2). 
For the electron-phonon term, taken from Eq. Al, the net wall energy due to both atoms is 

AE a = -a1R\\2 x 1 - + |1 + ( l -f - 0| - 4). 
Adding up the excess energies of these two cells, the wall energy is found to be 

E a w = atR 2 -(^)r R 4 

= (y^p ~^) Eo = 1A8E o(Ve P - 0.186). (B2) 

This procedure can be repeated for the thicker walls of Figs. 13b, 13c. In each case, the sum involves 
one additional cell, and there is now an internal degree of freedom, the tilt angle 9 of the cells within the 
wall. For the one-cell wall, Fig. 13b, 

„h E / . „ sin 4 8 + (1 + sin8) 4 sin 2 29 ,3nri er) . r _ „ „ , „ „ n > 

E b w = 2cos 4 9 - 3 + i- '— + — - — t] 2 + (—f*)[U - 2sin6 - \10cos 2 9 - 3 - 2sin6\]) 

4 V 8 2 8 

= 2.65£ (?7ep + 0.047?72 - 0.133), (S3) 

where the final numerical form is for the symmetric wall position, 6 = 45°. Actually, the term in rj ep is 
minimized at 9 = 0, and so dominates the wall energy that the total E h w also has a minimum at 9 — as 
long as ri e p > 1. In fact, this limiting case of 9 = corresponds exactly to the abrupt wall situation, Fig. 
13a. It is a check of the calculation that Eq. B3 agrees with Eq. B2 in this limit. 
The same calculation may be repeated for Fig. 13c. The wall energy is 

En / „ 0nn sin 4 9 + (sin9 + cos9) 4 + (1 + cos9) 4 9nn 3irri er) . n n 

E w = -^(- 3 - sin 2 29 + i ^ i '— + sm 2 2^ 2 + -^-[9 - 2cos6 - 5cos29}) 

= 2.80E (r] ep + 0.067?72 - 0.160), {BA) 

with the last line holding for 8 = 30°. No absolute value signs appear in Eq. B4, since it only makes sense 
for 9 < 45° (otherwise the wall tilt is no longer monotonic). Again, the wall energy actually minimizes at 
9 = 0, in which case it agrees with Eq. B2. Thus, the minimum wall energy corresponds to the abrupt 
wall, Fig. 13a. The energy, Eq. B2, should be compared to the continuum theory, Eq. 20. It is seen that 
the continuum theory underestimates the absolute energy by a factor of about two, but correctly predicts 
that the wall thickness is of the order of a single unit cell. It should be noted that the calculations of this 
Appendix employ the full, discontinuous form of the electron-tilt coupling (solid line in Fig. 14a) rather 
than the smooth approximation (dashed line). 

Note that the present model does not fully account for the complications of corner-sharing octahcdra. 
The distortion of, e.g., Fig. 13a is not consistent with a tilting of rigid octahedra. For instance, in Fig. 13a, 
the top atom to the left of the domain wall is tilting parallel to the wall (the upper apical oxygen is tipped in 
the positive y direction), so the in-plane oxygen in the right corner of the cell (positive x direction) would be 
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untilted, for a rigid octahedron. On the other hand, the apical oxygen of the top cell to the right of the wall 
is tipped in the negative x direction, which would require the same in-plane oxygen to tip below the plane. 
In the actual corner-shared configuration, neither tilt can be a pure x or y directed tilt. This distortion in 
turn will cause cells further from the domain wall to be tilted, and may lead to distortions of the shape of 
the octahedra. 

However, a more realistic model, allowing for these distortions, would be considerably more involved. 
For instance, it would include variations in the Cu-apical O distance, with the ensuing relative change 
of the Cu d 3z 2_ r 2 — d x 2_ y 2 level splitting. While these effects may play an important role in the detailed 
characterization of the domain wall and its excitations, it is hoped that the present simplified model captures 
the essential nature of the structural nonlincarities of the model. 











Table I: Parameter Values (from Ref. [13]) 


fit 




Vep 


V2 







7ioo 


Vl 




(meV) 












(m/s) 


10 


6.4 


5.6 


0.7 


68 


-3.5 


-0.22 


470 


11 


4.7 


9.6 


-2.0 


65 


-6.2 


-0.23 


400 
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Figure Captions 

Fig. 1 Proposed domain wall between two LTT domains, for (a) (110) wall or (b) (100) wall. Circles represent 
apical oxygens in layer above Cu02 plane. Arrow indicates direction of displacement of apical oxygen 
due to tilting of corresponding octahedron. Dashed lines indicate limits of domain wall, within which 
tilt direction is variable. Dotted lines indicate two columns of inequivalent atoms, A and B. 

Fig. 2 Phase portrait of Eq. 9 for fast waves (v > v\), with (3 > 0, 7 > 0. Lines with filled circles indicate 
separatrices. 

Fig. 3 Phase portrait of Eq. 9 for (3 < 0, 7 < 0. 

Fig. 4 Phase portrait of Eq. 9 for < 0, 7 > 0. 

Fig. 5 Temporal (or spatial) evolution of fast waves (corresponding to Fig. 2), assuming (3 = 7 = 0.2, a = 0, 
and yo — 2, 1, 0.82, and 0.8, in order of decreasing slope; dashed line is yo — 2 curve, with expanded 
horizontal scale (upper axis). 

Fig. 6 Oscillatory wave solutions, which can equally correspond to either fast waves (using left-hand axis) or 
slow waves (right-hand axis). Parameters (for slow waves) are (3 = 7 = —0.1, yo = 0, and a^/ir = 0.8, 
0.4, 0.2, 0.02, and 0.002, in order of increasing amplitude. For clarity, only the first two cycles of each 
mode are illustrated. 

Fig. 7 Orthon solution, assuming (3 = 7 = —0.1. (a) Phase of tilt; (b) cos2(j), proportional to the 'LTT strain', 

e_ = en — e 2 2; (c) sin2(j), proportional to the orthorhombic shear strain ei2- 
Fig. 8 Tetron solution, with same parameters as Fig. 7. 

Fig. 9 Family of tetron solutions, illustrating discontinuous jump in a for 7 > j c — 0.25. For all curves, 
[3 = -1. Values of 7 arc, from right to left, 0.05, 0.1, 0.2, 0.251, 0.26, 0.3, 0.4, 0.8, 2, and 10 (7 = 100 
is nearly indistinguishible from 7 = 10.) Filled circles show limits of jump in a. 

Fig. 10 Local energy per unit cell of waves of Fig. 6, assuming ij ep = 2, 772 = 0. For simplicity, only one half 
cycle of period (corresponding to a single kink) is illustrated. 

Fig. 11 Local energy per unit cell of single (100) soliton, using parameters appropriate to LSCO: r/ ep — 8, i] 2 = 0, 
Eq = b.hmeV . Solid line = total energy Eq. 17; dotdashed line = contribution of (3 e ; long dashed line 
= contribution of j e ; short dashed line = contribution of S e ; dotted line = contribution of — 3</>' 2 /8. 

Fig. 12 Integrated domain wall energy E w per column, as a function of wall spacing, d w , both for orthon (solid 
lines) and tetron configurations (dashed lines), for i] ep = 8 or 2. 
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Fig. 13 Microscopic domain wall configurations, for the calculations of Appendix II, for three different wall 
thicknesses. 

Fig. 14 (a) JT energy, H epn . Solid line = Eq. A2b, v = 0; dashed line — Fourier approximation, retaining only 
the n = 0, 1 terms, (b) dH epn /d<p n . Solid line = Eq. A4b, v = 0; dashed line = Fourier approximation, 
n = 1 term only; dotted line = two term (n = 1,2) Fourier representation. 
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